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Abstract. Canonical coupling between classical and quantum systems cannot result in 
reversible equations, rather it leads to irreversible master equations. Coupling of quantized 
non-relativistic matter to gravity is illustrated by a simplistic example. The heuristic derivation 
yields the theory of gravity-related decoherence proposed longtime ago by Penrose and the 
author. 



1. Introduction 

There is at least one classical dynamic system whose quantization is still problematic, that 
is gravity. Experimental evidences up to now have not indicated that gravitation should be 
quantized. Methods of quantizing the classical equations of gravity have failed to become 
generally accepted. So, gravitation could happen to be classical. Classical dynamics of 
gravitation would couple to quantum dynamics of other fields. This mathematical problem 
is not at all trivial. 

Hybrid quantum-classical systems are important not only for gravity but for molecular and 
nuclear systems as well where we might need to go beyond the mean-field approximation. 
Hybrid dynamics has obtained certain theoretical importance in foundations, in cosmology, in 
measurement problem. A very incomplete list of related works [Il[2l[3l[ll[5l[6l[71[8l[9l [TTl [12] 
shows the diversity of motivations. 

The particular goal of the present work is to re-derive a former heuristic gravity-related 
decoherence model. According to it, the Newton gravitational field (p possesses a universal 
uncertainty 5(j) which is modelled by a zero-mean white-noise with spatial correlation |13l [T3] 

Wr,t')<^0(^,t))„,_ = ^<J(t'-t), (1) 

where G is the Newton constant. This noise causes a certain universal decoherence between 
'macroscopically' different mass distributions. The following quantum master equation can be 
derived [HI [T5] for the density matrix pQ of the quantized system: 



where Hq is the Hamiltonian, Hq is the standard Newton pair-potential, /(r) is the operator 
of spatial mass density. (We have omitted the integral volume elements dr, ds and we keep this 
shorthand notation throughout our work.) The decoherence mechanism encoded in Eqs. ()ll2p 
coincides (upto a missing factor 1/2) with Penrose's proposal [T6l [T7] 

^ <y r^s I ' 5 1 

for the rate of decoherence between two different mass distributions / and /'. 

Sec. [2] introduces the mathematical model of hybrid dynamics which will be applied to the 
matter-gravity quantum-classical hybrid system in Sec. [3] including a re-derivation of Eqs. ()ll2p . 

2. Hybrid dynamics 

We are going to discuss a possible mathematical model for hybrid composite systems consisting 
of one quantum and one classical subsystem. 

2.1. Hybrid state, expectation values 

The generic state of a hybrid system is the hybrid density 

p = p{q,p) (4) 

which is a non-negative matrix in the Hilbert space of the quantum subsystem, depending on 
the phase space point {q,p) of the classical subsystem. Its quantum marginal 

PQ = j PiQ^P)dqdp (5) 
is the usual normalized density matrix of the quantum subsystem, while its classical marginal 

Pc{q,p) =ticp{q,p) (6) 

is the usual normalized Liouville density of the classical subsystem. (We mention the notion 
of conditional state PQ\qp = p{q-,p) / Pc{q-,p) of the quantum subsystem; conditional state of the 
classical subsystem wouldn't make sense in general.) 

The interpretation of the hybrid state follows from the interpretation of density matrices and 
Liouville densities, and goes like this. Consider an arbitrary phase-point-dependent Hermitian 
matrix A{q,p), i.e.: a hybrid observable. Its expectation value in the hybrid state p is defined 
by 

(^(9,P))/3 = tr j A{q,p)p{q,p)dqdp . (7) 

The expectation value has the usual statistical interpretation: if we measure the value of the 
hybrid observable A{q,p) repeatedly on a large ensemble of identically prepared hybrid systems 
characterized by the hybrid state p, the average value of the measured outcomes will tend to 
the above calculated expectation value. 

Once we have specified the general features of hybrid states, we can go on and learn about 
their dynamics. 



2.2. Dime, Poisson, Aleksandrov brackets 

Let us begin with two independent systems, one is quantum and the other is classical. Our 
quantum system's state is described by the density matrix pQ. The evolution is governed by the 
von Neumann equation with Hamiltonian Hq: 

^PQ = -^[Hq,pq] = (Hqpq - pqHq) . (8) 

For our classical system, we have chosen the classical canonical Liouville theory because of its 
best formal match with the quantum theory. The state of the classical system is described by 
the normalized phase space density pc = Pcil^p)- The evolution is governed by the Liouville 
equation with Hamilton function He = Hc{q,p). 



d_ ^ \ =sr^ ( dpc _ dpc dHc \ 



(9) 



where q = {qi,q2, ■ ■ ■ ,qn) and p = {pi,P2, ■ ■ ■ ,Pn) stand for the canonical pairs of coordinates 
and momenta. The classical structure { , }p is called Poisson bracket, the quantum structure 
— , ] is called Dirac bracket. 

Now we can take the above quantum and classical systems to construct a single composite 
hybrid system. We can easily form the hybrid composite state: 

p = p{q,p) = PQpc{q,p) , (10) 

as well as the hybrid Hamiltonian: 

H = H{q,p)^HQ + Hciq,p) . (11) 
We write the hybrid evolution equation in the form 

^^p=-^[H,p] + {H,p}p, (12) 

which is equivalent with the previous evolution equations ([5]) and Q for pq and pc separately. 

Of course, we would like to introduce interaction between the quantum and classical 
subsystems, i.e., we introduce the hybrid interaction Hamiltonian Hqc = HQc{q,p) which 
is a certain Hermitian matrix that also depends on the canonical coordinates of the classical 
subsystem. Now our total hybrid Hamiltonian reads 

H = H{q,p) =Hq + Hc{q,p) + HQc{q,p) , (13) 
and Aleksandrov ^ proposed the following evolution equation: 

j^p =-'-[H,p\ + Herm{F,p}p . (14) 

The structure on the r.h.s. is called the Aleksandrov bracket, being the simplest generalization 
of the naive interaction- free structure ()12p . 

The Aleksandrov equation ()14p does not preserve the factorized form (jlOp of the hybrid state 
p{q,p), which is the legitimate consequence of the interaction. But a serious drawback arises. 
As pointed out by Boucher and Traschen p], the positivity of p{q,p) may not be preserved, 
the Aleksandrov equation is fundamentally defective whereas it can be a useful approximation 
beyond the mean-field model. 



2.3. Blurring Dirac+Poisson 

However disappointing it is, the algebraic unification of classical and quantum dynamics runs 
into difficulties [2l[3l[ll[5]. Perhaps the deepest reason is that the notion of deterministic classical 
trajectories of one subsystem becomes lost under the influence of the other subsystem which is 
subject to quantum uncertainties. I proposed a possible remedy long time ago [3j, another 
approach was shown together with Gisin and Strunz [6]; further mathematical structures of 
hybrid dynamics appear from time to time [3 El O [TUl [TTl [T5]. A comparative analysis is 
missing. The present work follows the method of Ref . [3] . 

In the hybrid Hamiltonian (jl3|) . the interaction can always be written in the form 

Hqc = Y.Jq'Jc, (15) 

r 

where Ref. [3] called Jq and Jq as 'currents' that build up the hybrid coupling. This time, with 
one eye on the forthcoming application in Sec. O we write the hybrid interaction this way: 

^QC = E/V- (16) 

r 

Let us blur the interacting 'currents' by auxiliary classical noises 5f^ and 6(p'^: 

H^t" = E {r + ^fii)) i^" + '^<^'(*)) • (17) 

r 

Accordingly, we replace the hybrid Hamiltonian ()13p by the noisy one: 

and we replace the Aleksandrov equation (fT^ by its noise-averaged (blurred, coarse-grained) 
version: 

-^p= ,p] +Herm{i?™,p}p\ . (19) 

Ut \ n I noise 

Suppose the noises bf and b(\) are independent white-noises of zero mean, with the following 
auto-correlations : 

('^r(W(t)U,e = ^Q'^(*'-*)> (20) 
{b^\t')W{t))^^^^ = D}?5{t'-t). (21) 

The two correlation matrices Dq and D'^j can be chosen in such a way that the noise-averaged 
Aleksandrov equation (fTUj) guarantee the positivity of the coarse-grained hybrid density p. The 
price to pay is the loss of reversibility: the mathematically correct hybrid equation will be a 
certain irreversible master (kinetic) equation. 

2.4- Hybrid master equation 

By hand, we added a noisy part ()17p to the total hybrid Hamiltonian ()13p : 

fjnoise ^ ^ ^ Y.^f5(t/' + f6f) . (22) 

r 

Compared to Eqs. (|17|18p . the term 6f6(p has been omitted because it cancels from both the 
Dirac and Poisson brackets. In the evolution equation ()19p . the 1'^* order contribution of the 
noisy part is vanishing on average. In 2"*^ order, the term fd(f> adds the structure — [/, [/, p]] and 



the term (p6f adds the structure {(p, {(p,p}p}p to the r.h.s. of the naive Aleksandrov equation 



|p = -^[H, p] + Herm{i?, p}p - -L ^ D^Jif, [f , p]] + {'^^ {'/'^ pIpIp • (23) 

This hybrid master equation, i.e.: the noise-averaged Aleksandrov equation (fT9]l . preserves the 
positivity of p if the correlation matrices satisfy 

DqDc > ^ . (24) 

The proof was given in Ref . [3] (it has remained largely unnoticed) . 

The first irreversible term on the r.h.s. of the hybrid master equation ()23p is a typical 
decoherence term dephasing the superpositions between different configurations of the 'currents' 
/. The kernel Dq plays the role of decoherence matrix. The second irreversible term concerns 
the classical subsystem, yielding diffusion in its phase space. 

3. Gravity-related decoherence 

Hybrid dynamics will be applied to the interaction of quantized massive particles with classical 
gravitational field. 



3.1. Hybrid master equation for matter plus gravity 

The Hamiltonian Hq will stand for the quantized non-relativistic matter. We'll need its mass 
density field operator /, which should be the non-relativistic limit of Too/c^, the 00 component 
of the Einstein energy- momentum tensor divided by (c is the speed of light). For point-like 
particles of mass m and position operator x, the mass density would be 

/(r) = ^ m5{r — x) , (25) 

where summation extends for all massive particles. In fact, a certain short-length cutoff is always 
understood, for details of this issue see Ref. [TSl [THl [2D] and references therein. 

As to the classical gravitational field, we start with a scalar relativistic model and we take the 
non-relativistic limit c — )• oo afterwards. We define the Newton potential (f> = ^c^{goo ~ 1) where 
goo is the 00 component of the Einstein metric tensor. The canonical coordinates q of Sec. [2] will 
be played by the Newton field (/>(r), the conjugate canonical momenta p will be denoted by the 
field (,{r). The Hamilton function(al) can be chosen as 

^c(0,O = /(2vrGc2e^ + i^) • (26) 

To construct the hybrid interaction, we couple the quantized mass density / to the Newton field: 

HQci^) = I f{r)<P{r) . (27) 

Note, in parentheses, that Hqc does not depend on the canonical momenta Looking back 
to notations in Sec. [21 we recognize the identity of this interaction with (jl6p provided we 
replace J2r'f ''^'P'' by jr-, f{f)i(t^{f)-, respectively. Therefore the equations of Sec. [2] apply here 
conveniently. Let us form the total hybrid Hamiltonian (jl3p : 



(28) 



Write down the hybrid master equation ([23|) for the state p{(t>,i) of our hybrid system: 
j^P= -^[i?,/5]+Herm{i?,p}p 

For preserving the positivity of the state p, the decoherence and diffusion kernels should satisfy 
the condition (|24p . Below, we propose a heuristic resolution of the remaining ambiguity as to 
the concrete form of Dq and Dq- 

3.2. The gravity-related decoherence matrix 

In the Newtonian non-relativistic limit c — >■ cx), the mean- field Poisson equation is valid. In our 
notations ([7]): 

A (c/>(r))^ = 47rG (/>)), . (30) 

Let us recall the definitions ()20p and ()2ip of the correlation matrices (kernels) Dq and Dc, 
respectively. If we imposed the Poisson equation iS.5(j) = ^-nGSf on the fluctuations 5f, 54) in 
Eqs. ()20|2ip . respectively, then the above two correlation kernels would become related by 

AA'Z)c;(r, r') = {^7:GfDQ{r, r') . (31) 

(We are aware of a conflict with the derivation of the hybrid master equation where we considered 
5(j) and 5f uncorrelated. This is conceptual problem, it remains harmless for our resulting hybrid 
master equation but it may become relevant for later developments.) From the inequality ()24p . 
the minimum blurring corresponds to DqDq = h? /A. Hence the unique translation invariant 
solution reads 

Dc{r,s) = '^- (32) 

2 |r — s| 

This is the correlation ([1]) of universal gravitational fluctuations \\.2>\ [H] apart from the numeric 
factor 1/2. 

3.3. Reduced quantum master equation 

Let us invoke the expansion ^ of the Poisson bracket and calculate the full expansion of the 
r.h.s. of the hybrid master equation ()29p : 



|p=-i[fl«,pl-i/0M[/(r),pl 



We are going to derive the reduced dynamics of the quantized matter. Its density matrix is the 
marginal ([5]) of the hybrid state: 

PQ= j p{<t^,imVi. (34) 

Let us try to obtain a closed master equation for pq. Integrate both sides of (j33p by Vcp and 
V^. Most terms on the r.h.s. cancel and we are left with 

j^PQ = -^[Hq,Pq] -^j4>{r)[f\r),mm-^jDc{r,s)[f{r),[Hs),^ , (35) 



where /5((/)) = / /5((/), 4)2?^. Suppose the post-mean-field Ansatz: 



J Acf){r)p{(j))V(j) = 47rG Herm/(r)/5Q , (36) 

which is stronger than the standard mean- field equation ()30p . From this Ansatz, we express the 
functional integral on the r.h.s. of ([35]): 

J (t){r)p{(t))V(t) = -GHerm ^ y^^PQ ■ (37) 
Subsequently, the complete term in question becomes Hamiltonian: 

where Hq is the well-known Newtonian pair-potential: 

The Eq. (j35p has thus led to the following closed master equation for the evolution of the 
quantized matter: 

= -^[Hq + Hcpq] -^J^Dc(,r,s)[f{r), [/(s),/5q] , (40) 

which, invoking Dq from ()32p . is the master equation ^ \14:\ [15] apart from an additional 
numeric factor 1/2 (cf. also Ref. [21]) in front of the decoherence term . 



4. Concluding remarks 

We have constructed a possible consistent hybrid dynamics coupling classical gravity to 
quantized non-relativistic matter. We were able to re-derive the former heuristic equations 
of gravity-related universal decoherence. 

We admit a certain elusiveness of our procedure. In order to borrow self-dynamics to the 
Newtonian field (p, we attributed the status of relativistic scalar field to (p, with the non- 
relativistic limit to be taken on the hybrid master equation. On the other hand, the proposed 
hybrid master equation assumes the Markovian approximation since it contains classical white- 
noise fields. Coexistence of relativistic and white-noise fields may be an issue in itself. Therefore 
our hybrid master equation is heuristic rather than exact, and it remains so until, e.g., we 
construct some explicit solutions at least. The non-relativistic limit c ^ oo was treated formally, 
without check for its existence. This is why we keep talking about the post-mean-field Ansatz 
rather than the post-mean-field equation despite the heuristic proof presented in the Appendix. 

In any case, the formal re-derivation of the gravity-related decoherence master equation from 
a certain closed set of equations is reassuring and indicates the possible relevance of hybrid 
dynamics in foundations |22^ |23] . 



Acknowledgment s 

This work was supported by the Hungarian Scientific Research Fund under Grant No. 75129. 



Appendix 

We ought to show that the post-mean-field Ansatz ()36p is consistent with the non-relativistic 
limit c — )■ cxD of the hybrid master equation ()29p . We begin with the heuristic derivation of the 
standard mean-field equation ()30p . 

From the hybrid master equation ([33|) and the definition of expectation values ([7]), calculate 
the time-derivatives of the classical mean fields and {£,)p: 

±{^{r))^ = A7rGc\ar))p, (-1) 

Prom the first equation wee see that {S,)p — > in the non-relativistic limit. It also means that 
the l.h.s. of the second equation vanishes, hence we have: 

(Ac/>(r)-47rG/(r))^ = 0. (.3) 

This is just the standard mean-field equation (j30p . 



Let us upgrade it into the post-mean- field Ansatz ()36p . To this end, we calculate the time- 
derivative of the expectation value of the simple hybrid quantity A(j) where A itself is an arbitrary 
constant Hermitian matrix: 

= 47TGc'{Aar))p + 0{l) . (.4) 

Here 0{1) stands for terms that does not contain the diverging factor (P. In the non-relativistic 
limit, the quantity (A^(r))^ must vanish. Since A can be any Hermitian operator, the 'partial' 
expectation value, i.e.: the expression d?]) without the trace and A, does already vanish: 

I ^{r)p{^,OVm = 0. (.5) 



-(i0(s)F(O)p = 4^Gc2(iC(r)F(0)p + 0{1) , (.6) 



-(i^(r))^ = - (AAcPir) - 4^GHermi/(r))^ (.7) 



In general, the hybrid dynamics (|33p yields 

d 

m 

where -F(^) can be a generic functional of Repeating the previous arguments, we conclude that 
the 'partial' expectation value of ^{r)F[^) must vanish for c — >■ oo. Since F[^] is an arbitrary 
functional, we come to the conclusion that, in the non-relativistic limit, E,{r) is identically zero. 
Now we derive the time-derivative of {A^[r))p from ()33p : 

+ '^{\tiQ + jmf{s),M{r)\-\ jDc(n,,)\hu)Ah>).kir)\) ^ 

As we said, in the limit c ^ oo we can set ^ to zero identically for all time and position therefore 
the above equation reduces to: 

(iA0(r) - 47rGHermi/(r))p = . (.8) 

The l.h.s. vanishes for all A, hence we can remove the trace and A from the hybrid expectation 
value, yielding: 

'A(t){r) - 4^GHerm/(r)) i)V^Vi = . (.9) 



This is exactly the post-mean-field Ansatz (|36p if we invoke the definition ([5]) of pq. 
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